Freshwater Biology

The prevalence of nitrogen limitation and nitrogen-phosphorus co-limitation (henceforth referred to as nitrogen-related limitation) in Norwegian lakes and their relationships with atmospheric nitrogen deposition, climate, dissolved organic matter (DOM), and catchment characteristics were assessed across space and time. Routine monitoring data from 1,529 lakes in the national Vannmiljo database were analyzed for two multi-year periods (1995-2009 and 2010-2025). Limitation was inferred using the molar NO--N/TP ratio as an indicator of dissolved inorganic nitrogen availability. Nitrogen-related limitation was widespread in both periods and exhibited strong regional structure, with highest prevalence in northern regions and lowest prevalence in southwestern Norway. Overall prevalence increased from 31% to 38% between periods, with significant increases in western regions. Regional-scale models identified climate, forest cover, DOM, agriculture, and atmospheric nitrogen deposition as predictors of limitation probability, whereas study period per se and bog/peatland cover were not significant. At the local scale, atmospheric nitrogen deposition and DOM were the only consistent predictors, with substantially lower explanatory power than at the regional scale. These results indicate that large-scale environmental gradients play a major role in shaping nutrient stoichiometry in Norwegian lakes. Because the monitoring dataset primarily represents lakes affected by human activities, the findings are particularly relevant for water management. The widespread occurrence of nitrogen-related limitation suggests that nitrogen availability may influence phytoplankton growth in many systems and that dual-nutrient management strategies addressing both nitrogen and phosphorus may be required in many regions.

O_LIEnvironmental DNA (eDNA) metabarcoding is an emerging tool for biodiversity assessment in freshwater systems, offering high-resolution insights into community composition. Here, we apply eDNA metabarcoding to evaluate the ecological impacts of Eurasian beaver (Castor fiber) activity within a seminatural enclosure in the Scottish Highlands. C_LIO_LIWe collected seasonal water samples from nine sites, six influenced by beaver dams and three control sites with no evidence of beaver engineering, across a 40-hectare enclosure. Samples were analysed for vertebrate and macroinvertebrate diversity using established 12S and COI markers. C_LIO_LIVertebrate alpha diversity did not differ significantly between beaver and control sites, likely reflecting the small spatial scale and low species richness of upland Scottish streams. However, community composition differed significantly between treatments, especially for fish (PERMANOVA, R2 = 0.55, P < 0.001), with beaver-influenced sites dominated by three-spined stickleback and control sites by brown trout. Macroinvertebrate communities showed a 78% increase in gamma diversity in beaver-modified habitats relative to controls. Species composition varied strongly with beaver presence (PERMANOVA, R2 = 0.29, P < 0.001), likely due to the creation of lentic-lotic mosaics and associated microhabitat diversity. Seasonal variation was significant in both taxonomic groups, with the lowest species richness and highest community dispersion observed in summer, probably reflecting hydrological and temperature-driven dynamics in eDNA production and transport. C_LIO_LIOur findings reinforce previous evidence that beaver dam-building activity enhances beta diversity in headwater systems. Additionally, we demonstrate that eDNA metabarcoding is a sensitive method for detecting spatial patterns in freshwater biodiversity associated with these activities at scales ranging from tens to hundreds of meters. These approaches could inform future monitoring strategies aligned with landscape-scale beaver management and reintroductions. C_LI

Philippine freshwater ecosystems are considered one of the most diverse ecosystems harboring numerous fish species. However, in the Philippines, these ecosystems are threatened by invasive species that potentially disrupt ecological balance. In this study, we focused on the vermiculated sailfin catfish Pterygoplichthys disjunctivus, an invasive aquarium species reported in several Philippine aquatic ecosystems. Despite its documented spread, its potential range under a rapidly changing climate remains poorly understood for the country. Hence, in this study, we utilized the MaxEnt model to predict its near-current and future habitat suitability in the Philippines. Using 11 reported occurrences, our model showed high predictive accuracy (AUC = 0.882{+/-} .034, TSS = 0.7394 {+/-} 0.154, SEDI = 0.971 {+/-} 0.019). Across the current and future scenarios, slope was the primary contributor (78.7% - 81.3%), followed by BIO 10 or the mean temperature of the warmest quarter(18% - 27.8%), and flow accumulation (0% - 5.2%). However, for the SSP126 scenario, BIO10 is projected to triple by 2050 (18 - 48%). Current projections identify high-risk regions, particularly central Luzon (Laguna de Bay and Lake Taal), the Cagayan River Valley, and portions of eastern Mindanao (Agusan Marsh and Lake Mainit). Sankey transition analysis confirms a high habitat stability rate (>73%) for high-suitability pixels in both SSPs, indicating persistent invasion risk. Overall, our study provides a framework for invasive species management and contributes to the conservation of Philippine aquatic ecosystems.

Telomere length (TL) is increasingly used in ecology as a biomarker of individual quality and environmental stress, yet research on non-model species with complex life histories remains limited. Because TL varies among tissues and across ages in a species-specific manner, identifying non-lethal tissues that reliably reflect whole-organism telomere dynamics is essential for longitudinal telomere studies in the field. This study aimed to evaluate tissue-specific TL in Japanese eel (Anguilla japonica), an endangered catadromous fish. We first mapped the chromosomal distribution of telomeric sequences using fluorescent in situ hybridization (FISH), the first application of this method in this species. We then tested whether muscle and caudal fin, which can be sampled easily and non-lethally, can serve as suitable proxy tissues for TL measurements in wild individuals. Relative telomere length (RTL) was quantified by qPCR in blood, brain, caudal fin, gonads, heart, liver, and muscle. FISH analysis confirmed telomeric repeats at all chromosomal ends, with only weak interstitial signals on three chromosomal pairs unlikely to affect qPCR-based estimates. A generalized additive mixed model and Wilcoxons signed-rank tests revealed significant inter-tissue differences: RTL was shortest in the brain and muscle and longest in liver, blood and caudal fin. Muscle and caudal fin RTL were significantly correlated with RTL in many other tissues, supporting their use as proxy tissues for longitudinal TL monitoring, including responses to environmental variation. Both total length and age were tested as explanatory variables for RTL, and the model including total length showed a better fit than the age-based model. Non-linear relationships between RTL and total length observed in several tissues suggest physiological shifts associated with growth and sexual differentiation. Overall, these findings advance understanding of telomere dynamics in eels and establish muscle and caudal fin as suitable tissues for repeated, non-lethal TL assessment in ecological and conservation contexts.

Microbial communities are central to the biogeochemical cycling of nutrients, critically shaping ecosystem functioning and influencing climate change mitigation. Mangrove ecosystems are among the most important global carbon sinks that enable large amounts of carbon to be sequestered and stored. However, gaps persist in understanding the fundamental aspects of microbial-driven carbon cycling in these environments. This research explores the microbial taxonomic and functional diversity related to carbon cycling in selected tropical mangrove sediments across various locations and depths. Sequencing data analyses based on the 16S rRNA gene revealed distinct microbial community composition but conserved predicted functions across the different mangrove locations. Depth was a strong influence on the functional composition, with carbon-related pathways and metabolic strategies differing between top and bottom sediments. Putative functional gene abundance analyses revealed that carbon fixation processes were among the top carbon-related pathways, suggesting the key role of mangrove microbial communities in sustaining long-term carbon storage. Within these communities, Desulfobacterota appeared as a primary contributor to carbon fixation, while Chloroflexota played a significant role in carbon metabolism and methane cycling. Co-occurrence network analyses also revealed that these microbial groups were among the keystone taxa in mangrove sediments. Our study adds on to the body of knowledge on the mangrove microbiome and their carbon metabolic processes, which helps to improve strategies for managing and leveraging these vital carbon sinks.

AO_SCPLOWBSTRACTC_SCPLOWEnvironmental heterogeneity across freshwater systems often promotes phenotypic variation, yet disentangling environmentally induced variation from heritable differentiation remains a central goal in evolutionary ecology. We investigated the geographic distribution and morphological differentiation, and heritability of shell traits among populations of the freshwater lymnaeid snail Pectinidens diaphanus in Patagonia. Extensive field surveys across 196 freshwater sites revealed that the species occupies a broad range of lentic and lotic habitats and constitutes the only lymnaeid inhabiting southern Patagonia. While reproductive anatomical structures were conserved across populations, shell shape differed markedly among populations from contrasting habitat types, with population identity explaining nearly 50% of total shape variation. Populations from hydrologically unstable habitats (ponds and streams) exhibited more elongated shells and relatively smaller apertures, a pattern consistent with functional responses to hydroperiod variability and desiccation risk. To assess the heritability of this differentiation, we conducted a common-garden experiment across two generations. Shell shape differences between permanent- (lagoon) and temporary- (pond) habitat-derived populations persisted into the G2 generation reared under standardized laboratory conditions, indicating that the observed variation is not solely a response to local environmental conditions but includes a heritable component. Together, our findings demonstrate that P. diaphanus constitutes the sole lymnaeid across southern Patagonia, occupying a broader range than previously documented, and that populations show heritable shell differentiation potentially associated with contrasting freshwater habitats. By integrating large-scale biogeographic surveys with morphometric and experimental approaches, this study provides new insight into how habitat variation may contribute to ecological and evolutionary differentiation in freshwater gastropods.

Dormancy is a widespread adaptive strategy that allows organisms to survive in temporally varying habitats by suspending development and reproduction. Although environmental variability is expected to shape dormancy strategies, it remains unclear how differences in environmental variability and predictability influence both the production of dormant embryos and the termination of dormancy. We addressed these questions by comparing D. pulex and D. obtusa, two closely related species that inhabit environments differing in variability and predictability. We hypothesized that D. obtusa, which inhabits ephemeral environments, would exhibit a greater propensity for sexual reproduction and dormancy and would require stronger cues to break dormancy than D. pulex, which occurs in more permanent, predictable habitats. Consistent with our hypothesis, D. obtusa lineages produced significantly more males and ephippia than D. pulex when reared under identical laboratory conditions, indicating greater investment in sexual reproduction and dormancy. Contrary to our hypothesis, we found no difference in responsiveness to cues between the two species. Across species, embryos broke dormancy and hatched most readily after experiencing changes in cold and light, even if not experienced at the same time. In contrast, desiccation reduced the propensity to break dormancy. Together, these results indicate that species occupying more ephemeral environments invest more heavily in the production of dormant offspring, but that the environmental cues regulating dormancy termination appear broadly similar between species. This pattern suggests that while investment in dormancy may evolve in response to environmental variability, the mechanisms controlling dormancy termination are more conserved.

Early detection and monitoring of non-indigenous species (NIS) is crucial to prevent their establishment and to reduce ecological and economic impacts in coastal ecosystems. Traditional monitoring approaches, which rely largely on morphological identification of collected organisms, are often time-consuming and may fail to detect species that occur at low abundance, are morphologically cryptic, or are present in the form of inconspicuous life stages. DNA-based approaches, particularly those resorting to environmental DNA, have demonstrated high aptitude for biodiversity monitoring and biosecurity surveillance. By examining the genetic material from bulk community samples or released into the environment, DNA-based approaches enable the detection of species without the need for direct observation, thereby increasing detection sensitivity and expanding the scope of monitoring programs. Despite the rapid growth of its employment in marine monitoring, a global synthesis of the status and trends of DNA-based approaches for detecting NIS in this environment has been lacking. Here, we present such synthesis, based on 146 published studies employing DNA for NIS detections in coastal environments. Two main methodological approaches were used across the reviewed studies, namely DNA metabarcoding which was applied in 49% of studies, closely followed by targeted single-species PCR assays, used in 42% of the studies. A smaller proportion of studies (10%) combined both approaches, integrating broad community screening with targeted detection to improve surveillance efficiency. Globally, 752 NIS were detected across disparate taxonomic groups, with metazoans representing the largest proportion of detections (464 species), followed by Chromista (210 species) and Plantae (77 species). Among these, the most frequently detected taxonomic groups included Dinophyceae (Dinoflagellata), Teleostei (Chordata), Florideophyceae (Rodophyta), Polychaeta (Annelida), Copepoda and Malacostraca (Arthropoda), and Ascidiacea (Chordata). At the species level, several well-known marine invaders were recurrently reported, including Bugula neritina (Linnaeus, 1758), Styela plicata (Lesueur, 1823), Acartia (Acanthacartia) tonsa Dana, 1849-1852, and Botryllus schlosseri (Pallas, 1766), highlighting the ability of DNA approaches to detect widespread and established invaders across different regions. The mitochondrial cytochrome c oxidase subunit I (COI) gene was the most widely used genetic marker, reflecting its broad taxonomic coverage and extensive representation in reference databases, particularly for targeting Metazoa. Ribosomal RNA genes, particularly 18S and 16S rRNA gene markers, were also frequently employed to target a wider range of eukaryotic taxa. Regarding sampled substrates, water was by far the most analyzed substrate, followed by zooplankton and biofouling communities collected from man-made structures. Notably, approximately 31% of all NIS detections reported in the reviewed studies constituted new regional records. These results highlight the potential of eDNA for coastal monitoring but also underline important limitations. Persistent geographical, taxonomic, and methodological biases can affect detection outcomes, and reliance on single sample types or markers may increase false negatives - particularly critical for NIS early detection. Therefore, multi-marker and multi-substrate approaches are essential to improve detection reliability and support effective biosecurity strategies. As reference databases continue to expand and methodological protocols become increasingly standardized, DNA-based monitoring is likely to play a central role in future management and surveillance of biological invasions in coastal ecosystems. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=133 SRC="FIGDIR/small/722998v1_ufig1.gif" ALT="Figure 1"> View larger version (75K): org.highwire.dtl.DTLVardef@17948b1org.highwire.dtl.DTLVardef@193832dorg.highwire.dtl.DTLVardef@189033dorg.highwire.dtl.DTLVardef@33cddf_HPS_FORMAT_FIGEXP M_FIG C_FIG

O_LIThe frequency and intensity of heat events is increasing across marine and terrestrial ecosystems. Within the same ecological community, the relative exposure and sensitivity to heat stress may vary considerably among interacting species, like predators and prey. This can be especially true for species that interact at the aquatic-terrestrial interface, as well as for interactions between primarily nocturnal and diurnal species, making it difficult to predict how such communities will respond to habitat warming. C_LIO_LIThermal limit metrics such as CTmax are often assumed to equate with ecological death because such temperatures impair behavioral activity and/or physiological functioning. Prey that are diurnally active can be more frequently exposed to temperatures that approach CTmax compared to their nocturnal predators, which may use thermal refuges during the day. Yet the impacts of daytime heat exposure on nighttime predation risk remain unknown. C_LIO_LIHere, we compared the thermal environment, performance, and heat tolerance between the predatory blue crab, Callinectus sapidus and one of its prey species, the mangrove periwinkle Littoraria anguilifera in a tropical mangrove ecosystem. We examined how exposing prey to heat stress at and below their CTmax affected their capacity to avoid predation in the field at night when predation risk is highest. C_LIO_LIWe found that acute exposure to temperatures near CTmax during the day increased the prey species susceptibility to predation during recovery at night. Although both interacting predator and prey have high thermal tolerance, prey are exposed to conditions that already reach CTmax, suggesting that current extremes in temperatures may already be influencing vulnerability to predation in this ecosystem. C_LIO_LIOur results suggest that differential exposure to sublethal heat stress in diurnal prey relative to their predator, along with the subsequent impact of these exposures on predation risk, will play a role in shaping these interacting as climate warms. C_LI

Understanding risks of biological invasions associated with ballast water (BW) requires full understanding of the biodiversity transported in ballast tanks. Here we characterize the remarkable level of diversity that can be carried in the BW of a single vessel. To maximize our ability to capture BW diversity we: 1) utilized DNA-based methods to describe biodiversity, including both native and non-native taxa; 2) exploited multiple primer sets targeting multiple genomic loci with different expectations for taxonomic coverage; 3) assessed multiple tanks on a single vessel to capture different communities present in different tanks; and 4) sampled those tanks with far higher-than-usual replication both to improve representation of the diversity present and to enable statistical estimation of total richness. Using this approach, we found extraordinarily high diversity associated with a single vessel. Across all loci, we estimate a total of 272 taxa that can be assigned species names; looking more broadly at unnamed molecular operational taxonomic units, our estimates are between 425 and 742 individual taxa, depending on the locus. We confirm that only a fraction of this diversity would be captured with typical sampling efforts. We found that different loci capture different snapshots of biodiversity and that our ability to detect taxa of interest (e.g., non-native species) depends on sampling effort and genomic locus. Our results expand upon previous studies describing highly diverse BW communities and add to a growing literature that demonstrates the value of molecular methods for characterizing those communities and assessing the associated risk of non-native species introduction.

The maintenance of endosymbiosis in cnidarians depends on the tight regulation of host immunity, cell cycle, and nutrient exchange, yet how these processes are impacted by interacting environmental stressors remains largely unknown. To address this, we employed physiological metrics, gene expression analysis, microbiome characterization, imaging (NF-{kappa}B localization, endoplasmic reticulum ultrastructure, EdU labeling), and stable isotope tracing in the model sea anemone Exaiptasia diaphana to examine the effects of heat and nitrate on these regulatory processes, individually and in combination. Heat treatment led to NF-{kappa}B activation, proteostatic stress, suppression of nutrient exchange, decreased cell-cycle progression, and microbiome restructuring, with all effects more pronounced in symbiotic than aposymbiotic anemones. In symbiotic anemones, nitrate partially offset these heat-induced responses through sustained carbon translocation, suggesting that the presence of symbionts, in conjunction with elevated nitrate, can temporarily buffer host thermal stress. However, prolonged combined exposure resulted in holobiont failure. These findings reveal that while nitrate enrichment can transiently delay the onset of bleaching, it does not preserve the regulatory networks required for symbiotic stability -- underscoring the vulnerability of cnidarian holobionts to the compounding effects of warming and nitrate pollution.

Microplastics (MP) are widespread in aquatic ecosystems and pose a threat to freshwater biodiversity. While numerous studies examine physiological effects on aquatic organisms, less is known about how MP alter chemically mediated interactions that regulate predator-prey dynamics. Predator-induced defenses in Daphnia depend on detecting kairomones and represent an important form of adaptive phenotypic plasticity. Whether MP interfere with these responses, and through which mechanisms, remains unclear. Here, we show that polystyrene MP impair predator-induced defenses across Daphnia species by disrupting predator-cue-mediated plasticity at the behavioral, morphological, and molecular levels. In D. longicephala, chronic exposure to PS fragments weakened Notonecta-induced morphological defenses, whereas additive-containing PS fragments nearly suppressed defense formation and reduced body size. Consistent with these phenotypic effects, proteomic analyses revealed alterations in pathways related to molting and chitin metabolism, linking MP exposure to impaired defense formation. In D. magna, PS particles attenuated fish kairomone-induced diel vertical migration, with stronger effects for larger particles, consistent with reduced effective availability or perception of predator cues. Natural limestone particles caused only minor effects, indicating particle-specific rather than general particle-driven responses. Our findings demonstrate that MP can disrupt adaptive predator-prey interactions with potential cascading consequences for freshwater food webs.

Global changes in land use and nutrient cycling are transforming ecosystems at unprecedented rates, with significant consequences for infectious disease dynamics. Aquatic environments are particularly vulnerable because the interplay of habitat modification, nutrient enrichment, and biodiversity loss can drive pronounced changes in the community composition of food webs, including hosts and parasites. Yet, despite well-documented effects of habitat modification on aquatic communities and food webs, the mechanisms through which these changes influence infectious disease dynamics remain poorly resolved. This gap arises, in part, because it remains challenging to disentangle how multiple stressors interact to shape disease outcomes and quantify parasite levels and host densities from field-collected samples. Here, we illustrate two tools that might help address these challenges. First, highly sensitive droplet digital PCR can quantify infection loads even when the signal:noise ratio is low. Second, stepwise Bayesian path analyses can identify the direct and indirect pathways connecting land-use changes to infectious disease dynamics. As a case study, we examined cyclopoid copepods and their helminth parasite, Schistocephalus solidus, across 47 freshwater lakes on Vancouver Island, a region strongly shaped by commercial logging, including widespread clear-cutting of old-growth forests. Our results reveal a positive correlation between copepod density and deforestation, potentially mediated by associated changes in water quality and calanoid copepods, key competitors of the focal host. ddPCR enabled sensitive detection of extremely low parasite signals in field-collected copepods. We detected positive infections in only 19.5% of the lakes surveyed, highlighting the difficulty of assessing disease dynamics in natural populations. Nonetheless, this study highlights the challenges of linking land-use change to disease outcomes, while also demonstrating that sensitive molecular and statistical tools offer new ways to reveal these hidden connections.

The sunflower sea star (Pycnopodia helianthoides) suffered a catastrophic population decline across its range from 2013 to 2017 due to the devastating Vibrio pectenicida FHCF-3 driven sea star wasting disease (SSWD) pandemic with minimal signs of population recovery. The functional extinction of this apex predator across substantial parts of its range has created a need to identify and track the remaining intact populations. Environmental DNA (eDNA) approaches provide a simple, cost-effective, and non-destructive method for monitoring occurrences, and in some cases abundances, of marine species, consistently outperforming visual occurrence monitoring efforts in sensitivity, speed, and cost. Here, we designed, developed, and validated a P. helianthoides-specific eDNA assay to identify refugia, using both quantitative and digital droplet PCR approaches. We first generated the most comprehensive sea star mitochondrial genome reference database to date (n=93 taxa, n= 15 novel). We then used unikseq and Geneious bioinformatics software to identify the unique nad5 gene region and design a highly specific hydrolysis probe-based PCR assay. We validated the performance of this assay through laboratory, mesocosm, and field testing, demonstrating a highly specific and sensitive assay. In a field application of the new assay across regions in British Columbia, Canada, we found a positive correlation between P. helianthoides eDNA concentrations and biomass density, especially when appropriately accounting for spatiotemporal integration scales (R2=0.67). The eDNA assay provides a rapid and scalable tool for monitoring the sunflower sea star which has been proposed for listing as threatened under the U.S. Endangered Species Act of 1973. Molecular tools like the one presented here enhance management and recovery efforts not only by identification and monitoring of remnant wild populations, but also by helping to assess population level response and recovery following reintroduction efforts.

Aimto characterize the evolution of climatic niches during the diversification of the Phyllotis darwini species group, in order to assess the extent to which divergences involved in radiation were associated with patterns of conservatism or divergence of climatic niches, and whether the differentiation found among climatic niches correlated with species phylogenetic relationships. Locationsouth-central Andes, surrounding lowlands, and Patagonia, South America. Methodsspecies climatic niches were characterized by sampling contemporaneous precipitation and temperature conditions across occurrence locations and entire distributional ranges. Climatic niches were analyzed and modeled using multivariate statistics (PCA, PERMANOVA), a maximum entropy-based algorithm, and novel methods developed to explore levels of differentiation (niche overlap test) and divergence (niche divergence test) between realized and fundamental niches. Comparative phylogenetic methods were applied using a time-calibrated phylogeny and integrating climate niche data to estimate ancestral environmental niches within geographic and environmental spaces. Resultscomparisons revealed low levels of climatic niche overlap, both among species realized niches and among their fundamental niches, suggesting high levels of niche differentiation during the diversification of Phyllotis species. Quantifications of niche overlap further showed that observed differences among species lay primarily in the multidimensional nature of climatic niches, as unidimensional quantifications exhibited higher levels of overlap. Evolved differences among species climatic niches were better fitted to a Brownian motion model of evolution, but lacked phylogenetic signal and showed no significant association with species phylogenetic distances. Main conclusionslow levels of differentiation between ancestral climatic niches suggest that the early radiation of species in the Phyllotis darwini species group was promoted by geographic isolation, whereas the more recent diversification of extant species was accompanied by climatic niche differentiation, possibly involving local adaptation to regional ecoclimatic changes associated with Quaternary glacial cycles. The spatial separation of sister species, the complete divergence of their climatic niches, and the lack of phylogenetic signal in niche differences suggest a scenario of diversification in which divergences were prompted by the spatial isolation, but also by the divergent selection exerted by regional climatic differences.

The mechanisms driving host-symbiont associations across space and time in contemporary mutualisms can give insight into the capacity for symbiotic organisms to respond to environmental change. High specificity between partners can increase cooperation and facilitate efficient holobiont selection, whereas low specificity could reduce host benefit, but facilitate adaptive associations across heterogeneous environments. The present study explores specificity in natural populations of a cnidarian-algal model, Exaiptasia diaphana, across a latitudinal gradient to understand the genetic and environmental effects driving host-symbiont associations, and their relation to heritable and/or environmental symbiont acquisition. We found that symbiotic associations were extremely flexible in E. diaphana, regardless of transmission mode. E. diaphana were capable of associating with diverse symbiont communities across genetically identical hosts seeded with vertically transmitted symbionts, as well as across highly connected host populations which acquire symbionts horizontally. Host population connectivity was complex and unrelated to geographic distance, whereas symbiont community composition tracked the thermal gradient, potentially due to context dependent biotic interactions. These results indicate that in a flexible symbiosis, symbiont communities are environmentally-determined, suggesting the future of this symbiosis will likely depend on climate adaptation of symbionts.

Environmental DNA (eDNA) surveys are increasingly used to assess marine biodiversity and inform deep-sea environmental decision-making, including mineral resource management and fisheries oversight. Yet standard low-volume protocols inherited from coastal work may be inadequate at depth, and no quantitative framework links depth and ecosystem context to defensible filtration volume targets. We compiled 841 eDNA samples from eight expeditions across the North Atlantic, Wider Caribbean, and Pacific (surface to 4000 m) to quantify how recoverable eDNA scales with depth and surface productivity, and to derive depth- and productivity-aware sampling targets. Total eDNA concentration declined with depth as a power law, with attenuation exponents (b) modulated by surface productivity: most gradual in eutrophic waters (b = 0.67), intermediate in mesotrophic (b = 0.90), and steepest in oligotrophic systems (b = 1.25); volume-weighted models explained 66-88% of the variance. At a fixed extract-concentration target, required filtration volumes diverged ~7-fold between oligotrophic and eutrophic systems at 200 m and ~38-fold at 4000 m. Conventional Niskin sampling, therefore, undersamples deep-sea biodiversity, particularly in mid- to low-productivity systems. Among laboratory parameters, the assay-specific extract-concentration target exerted greater leverage on required volume than extraction efficiency or elution volume. Volume-aware sampling paired with optimized recovery should be routine in deep-sea eDNA surveys.

Guinea worm disease eradication efforts may benefit from environmental surveillance methods capable of detecting infected copepod intermediate hosts in aquatic habitats. We developed a three-dimensional, spatially explicit agent-based model to examine how ecological processes influence detection probability for a hypothetical water sampling method. The results show that surveillance sensitivity is shaped by the combined effects of larval diffusion, copepod density, and pond size, with interactions among these factors producing nonlinear relationships. Detection, in our model, was concentrated within a relatively restricted period after larvae matured to the infective stage and before dispersal and mortality reduced presence, indicating a limited spatiotemporal window for effective sampling. Surveillance performance peaked under intermediate dispersal regimes that generated sufficient spatial overlap between larvae and intermediate hosts, while both limited dispersal and excessive diffusion reduced detection by constraining encounters or diluting larval concentrations. Increasing habitat size reduced detection by diluting larval concentrations, but the magnitude of this effect depended on copepod density and dispersal dynamics, producing nonlinear and threshold responses rather than simple scaling with pond volume. Spatial and temporal patterns of detection shifted as larvae dispersed, with the most favorable detection periods occurring when both larval abundance and intermediate host encounters were elevated. These findings indicate that surveillance can be guided by local ecological conditions. When the timing of larval introduction is uncertain, effective surveillance requires repeated sampling over time to capture transient windows of detectability and the sampling will be less effective in very stagnant and highly mixed waterbodies. Overall, this study demonstrates how mechanistic modeling can support the design and interpretation of environmental surveillance strategies for Guinea worm eradication programs. Author summaryGuinea worm disease is close to eradication but confirming that transmission has fully stopped remains difficult because detecting infectious larvae in water is challenging. Transmission depends on freshwater copepods that become infected after ingesting Guinea worm larvae. These copepods are short-lived and unevenly distributed within ponds, and infected individuals may die before larvae reach the infective stage. As a result, environmental detection is inherently uncertain. We developed a three-dimensional agent-based model to simulate larval dispersal, copepod infection, and water sampling in a pond environment. The model shows that detection is constrained to a brief period when mature larvae and copepods overlap in space and time, and that this window depends strongly on local ecological conditions such as larval dispersal, copepod density, and pond size. Because infected copepods can be present outside these narrow detection windows, negative water samples do not necessarily indicate absence of transmission, highlighting the need for repeated, spatially targeted surveillance during the final stages of eradication.

Ocean warming is altering abiotic environments and biotic interactions experienced by marine organisms, where sensitive early developmental windows occur in biologically complex seawater communities. The impact of these interactions on developmental processes and fitness in hosts is not well understood, but likely contingent on the establishment of a host-associated microbiome. Here, we hypothesize that temperature and microbial exposure during embryogenesis influence larval microbiome assembly and host morphology. Strongylocentrotus purpuratus embryos were raised in low microbial richness (LMR) or high microbial richness (HMR) seawater at ambient (14 {degrees}C) or elevated (18 {degrees}C) temperature, then collected at 2, 4, and 6 days post-fertilization (dpf) following multiple feedings. Higher microbial diversity was observed in larvae that developed in HMR seawater when compared to LMR. Differences in relative abundances of dominant microbial families between seawater and larvae suggest some degree of host selectivity in microbiome assembly. Temperature did not strongly alter microbiome composition, but both temperature and microbial condition led to differences in larval morphology by 6 dpf, potentially due to enrichment of microbes with chemoheterotrophic functions. By linking how temperature and microbial communities interact with host development, we contribute novel insights into how early-life environmental conditions impact holobiont formation and morphology. One sentence summaryEarly developmental temperature and microbial conditions shape larval microbiome establishment and morphology.

Near-shore coral reefs in southern Guam (Mariana Islands) experience severe sedimentation, in particular during the wet season when rainfall and erosion are high. We sampled fragments of the reef-forming coral Porites lobata from opposite ends of a sedimentation gradient in Fouha Bay, southern Guam, during dry and wet seasons. Using DNA metabarcoding, we characterized the diversity and composition of P. lobata-associated Symbiodiniaceae and bacterial microbiome communities. As in many species of Porites, Symbiodiniaceae communities of P. lobata were dominated by variants of Cladocopium C15 with sites showing differences in Symbiodiniaceae communities attributable to variation in these Cladocopium C15 variants. Bacterial microbiomes of P. lobata were dominated by Endozoicomonadaceae, a family of putative coral bacterial endosymbionts involved in nutrient cycling. Site and seasonal differences in bacterial diversity and community composition were apparent. In close proximity to the mouth of the river draining into Fouha Bay, bacterial diversity was highest during the wet season when sedimentation is generally severe. Microbiome reorganization in response to sedimentation may explain this result, but we also found overrepresentation of bacteria associated with terrestrial origin close to the river mouth and/or during the wet season. Together these patterns highlight that coral Symbiodiniaceae and bacterial communities are both spatially and temporally structured in this disturbed system. IMPORTANCEThis study provides a time series dataset of coral-associated microorganisms, including dinoflagellate algae and bacteria, from a tropical bay impacted by sedimentation that results from upstream erosion of disturbed soils. Characterizing temporal patterns of coral-associated microbes provides insights into the dynamic nature of these communities. While microbiome variability across sites and seasons may be a result of acclimatization to different environmental conditions, we identified bacterial groups of putative terrestrial origin in sampled coral microbiomes that may have been exported from eroded soils to the near-shore reef. Considering that disturbed soils act as hotspots for the proliferation of potentially harmful substances, such as antimicrobial resistance genes, understanding microbial community connections at the marine-freshwater-terrestrial interface is an important step toward evaluating environmental impacts across connected ecosystems from ridge to reef.